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The nonlinear propagation of a circularly polarized electromagnetic (CPEM) wave 
in a strongly magnetized electron-positron-ion plasma is investigated. Two coupled 
equations describing the interaction between a high-frequency CPEM wave and the 
low- frequency electrostatic wake field are derived. It is found that the generation of 
the wake fields depends on the presence of the ion species. The wake field generation 
in turn leads to de-acceleration and frequency down conversion of the electromagnetic 
pulse. 
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I. INTRODUCTION 



Electron-positron (pair) plasmas are composed of electrons and positrons which have the 
same mass and opposite charge. Such pair plasmas appeared in the early Universe , and 
are frequently encountered in bipolar outflows fiets) in active-galactic nuclei 3^, [J, in micro- 

nnnr □ 

quasars [5i], in pulsar magnetospheres [g, l7|, l8|, l9j , in magnetars [10|, in cosmological gamma 
ray fireballs jll], in solar flares [l^, and at the center of our galaxy [l^. Multiterawatt 
and petawatt laser pulses interacting with solid density matters can create pair plasmas 
Q, Q, Q, Colliding electromagnetic pulses also generate pairs from vacuum 18 1. 

The pair plasmas at the surface of fast rotating neutron stars and magnetars are held in 
strong magnetic fields, while super-strong magnetic fields can be created in intense laser- 
plasma interaction experiments. Accordingly, the understanding of collective phenomena in 
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strongly magnetized pair plasmas has been a topic of significant interest 
Specifically, it is to be noted that in a magnetized pair plasma, we have new wave modes 
whose counterparts do not exist in an electron- ion magnetoplasma 23|]. Furthermore, both 
astrophysical and laboratory plasmas contain a fraction of ions besides the electrons and 
positrons. In such environments, the linear and nonlinear wave propagation characteristics 
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are modified 
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In this paper, we consider the nonlinear propagation of magnetic field-aligned circularly 
polarized electromagnetic (CPEM) waves in a pair-ion magnetoplasma, to show the possi- 
bility of exciting electrostatic plasma wakes by the ponderomotive force of the CPEM waves. 
Our manuscript is organized in the following fashion. In Sec. II we discuss the linear dis- 
persion relation, the group velocity and the group dispersion for the magnetic field-aligned 
CPEM waves in a pair-ion magnetoplasma. Section III considers the nonlinear interactions 
between finite amplitude CPEM waves and electrostatic plasma oscillations (EPOs), and 
presents the coupled equations. Nonlinearities associated with the electron and positron 
mass variations, as well as with the CPEM ponderomotive force are incorporated in the dy- 
namical equations for the modulated CPEM waves and the driven EPOs. Section IV presents 
the conservation laws for nonlinearly coupled CPEM wave EPO systems, and discusses the 
excitation of wakefields . Section V contains a brief summary of our investigation. 
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II. THE LINEAR DISPERSION RELATION 



We consider a magnetized quasi-neutral electron-positron-ion plasma, i.e. noe = noi + UQp, 
where noe is the background electron number density, nop is the background positron density, 
UQi is the background ion density, and the background magnetic field is given by Bq = Bqz. 
For a right hand circularly polarized electromagnetic (RCPEM) wave propaga ting along z 
with frequency to ^ tOpi = {A-KriQie^ /rriiYf'^ , we have the dispersion relation 
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where Upj = (47rnoje^/me)^/^, with j = e,p, is the plasma frequency, Uce = —eBo/nie = —oJcp 
is the electron gyrofrequency, e is the magnitude of the electron charge, and rUe is the electron 
mass (for left hand polarization uo + ujc'is replaced by uj 



IjJr 



35|). In the limit of a strong 



background magnetic field, i.e. oj ^ IcJcel, we can simplify ([T]) to obtain 
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where we have introduced the notation f2j = (47rnoje^/me)^/^ and used the background 
charge neutrality condition. From Eq. ([21) we obtain the group velocity 
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and the group velocity dispersion v'g = dvg/dk = {vg/k){l — Vg/c^). In the limit kc ^ Up/ucp, 



we obtain 



f p ~ c 
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and v'g ~ Vlf / 4:Uj1pk'^c, and we have neglected the root with negative sign of the group velocity. 
Thus, in this limit the RCPEM waves propagate with a speed close to the speed of light. 



III. NONLINEAR EVOLUTION EQUATIONS 



Next, we consider a modulated nonlinear CPEM wave (cu, fc) , where the ponderomotive 
force induces slowly {dt ^ oj) varying electrostatic oscillations. We still assume that the 
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motion is fast enough for the ions to be immobile. The evolution equation for the vector 
potential amplitude A in the Coulomb gauge can then be written as [36;] 



t{dt + Vgd,)A + lv'd',A- 



e,p 



rii 



e^uj^\A\ 



A = 0, (5) 



where the group velocity and dispersion are given above, and the indices e,p have been 
omitted for notational convenience. The electrostatic perturbations are described by the 
density fluctuation ni, the longitudinal velocity Vz and the electrostatic potential 0. These 
variables satisfy the continuity equation 



dn\ ^ dvz 
dt ° dz 



and the momentum equation 
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-e for electrons and g = e for positrons, together with 
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Here qF denotes the ponderomotive force, which can be written as 
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in a magnetized plasma. We assume that the temperature T is the same for electrons and 
positrons. 

Next, we consider a high-frequency wave propagating with the group velocity, such that 
the driven low-frequency perturbations to first order are stationary functions of the comoving 
coordinate z — Vgt. Together with the assumption uj <^ Uc, the ponderomotive force (Q then 
reduces to 
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where we have used VtflujuJc ^ 1 in the last step. Next we expand Eq. ([5]) in the small 
parameter uj/uJc to obtain 
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A = 0. 



As we can see, the correction term proportional to lojiOf. must be kept above, since there 
will be an approximate cancellation of terms proportional to ni/riQ and Correction 
terms proportional to the small parameter 1 — Vg/c that have been omitted in Eq. ffTTj) in 
general give an additional contribution to Eq. (|TT1) . However, unless we consider the regime 
where uj < ujp (in which case ion-motion must be included in the low-frequency dynamics), 
such terms are of higher order, and hence omitted. Next combining Eqs. ([6])- ([8]) we can 
express the density perturbations in terms of the ponderomotive force and the potential. 
Substituting this expression in Eq. ([5]), using oj <^ ojc and = T/mf, <^ (? to simplify the 
equation, we find 
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Next, combining Eqs. 



to solve for the scalar potential, we have 
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where uj'^^^ = cu^g + u^p. Integrating (fT3!) . dropping the term proportional to (f^ + Vp)/c^, 
and using the background charge balance we find 
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Qf2 + ^Ptot 
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The coupled equations ( fT2l) and ( fT4l) describe the nonlinear self-interactions of a high- 
frequency CPEM pulse combined with the generation of a low-frequency electrostatic wake 
field in a strongly magnetized electron-positron ion plasma. We note that there is no wake 
field generation in the absence of the ions. 
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IV. CONSERVATION LAWS 

In a coordinate system moving with the group velocity and with suitable normalizations 
(i.e. ^ = ujptot{z — Vgt)/c, etc.), the above system can be put into the generic form 



.dA d'^A 
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where a is a constant that determines the relative importance of the self-nonlinearity, in our 
case given by 
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Various aspects of coupled system of this type have been studied by Refs. 
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Following the presentation in 
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371] we will refer to the low-frequency field $ as the wake-field. 



The equation system (ITSl) and (fT6l) can be derived from a variational principle. Introducing 
the Lagrangian density 



^ 1 f ,,dA dA* 
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where the action functional is A{^, A, A*) = J CdC,dT, we obtain Eqs. (fTSll and ( fT6l) by 
varying $ and A* and minimizing the action as usual. The system posses the following 
three conservation laws: 



d 
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where W<s, = (1/2) [<I>^ + (9$/9^)^] . Within the framework of the WKB-approximation, 
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A\ d^ is proportional to the number of high-frequency quanta N, and thus Eq. f[T9|) 
describes the conservation of N. Using a quantum mechanical analogue, the process of wake 
field generation can thus be viewed as a parametric process where a high-frequency quanta 
with frequency u decays into a low frequency quanta Q and another high-frequency quanta 
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with frequency 00 — Q, conserving the total number of high-frequency quanta. Furthermore, 
Eq. fl20|) is the conservation equation for the Hamihonian. Although fl2Tl) does not look like a 
conservation law, it actually describes the conservation of energy. To see this, it is convenient 
to choose the positions ^± before and after the pulse passage, respectively. For this choice 
we note that the left-hand side of (12T]) is proportional to d/dr K\Ai^\^dK/ IAk^c/k = 
dAk/dr = {l/vg)dAuj/dT. Furthermore, since N is conserved, dAu/dr is proportional to 
the high-frequency energy decay rate. Finally, the difference in before and after the 
pulse passage is proportional to the energy transfer rate to the wake field, and thus we 
deduce that Eq. (l2Ti) implies energy conservation. 

V. SUMMARY AND CONCLUSIONS 

To summarize, we have considered the nonlinear interactions between magnetic field- 
aligned CPEM waves and electrostatic plasma oscillations (EPOs) in a pair-ion magneto- 
plasma. It is found that the EPOs are generated by the ponderomotive force of the CPEM 
waves only when the ions in a pair plasma are present. The present nonlinear wave- wave 
interactions provide the possibility of CPEM pulse localization and the generation of wake- 
fields. Localized EM pulses, in association with wakefields, can be identified in observations 
from astrophysical settings as well as from next generation intense laser-solid density plasma 
experiments where pairs and ions appear simultaneously. 

Acknowledgment s 

This research was partially supported by the Swedish Research Council (VR). 



8 



[a) ] Permanent address: Institut fiir Theoretische Physik IV, Fakultat fiir Physik und Astronomie, 

Ruhr-Universitat Bochum, D-44780 Bochum Germany. Also at School of Physics, University 
of KwaZulu-Natal, 4000 Durban, South Africa; SUPA Department of Physics, University 
of Strathclyde, Glasgow, Scotland; GoLP/Instituto de Plasmas e Fusao Nuclear, Institute 
Superior Tecnico, 1049-001 Lisboa, Portugal; Department of Physics, COMSATS Institute of 
Information Technology, Islamabad, Pakistan. 

[b) ] Also at Department of Physics, Linkoping University, SE-58183 Linkoping, Sweden. 

[1] W. Misner, K. S. Thorne, and J. A. Wheeler, Gravitation (Freeman, San Francisco, 1973). 
[2] M. J. Rees, The Very Early Universe, eds. G. W. Gibbons, S. W. Hawking, and S. Siklas 

(Cambridge University Press, Cambridge, 1983). 
[3] M. C. Begelman, R. D. Blandford, and M. J. Rees, Rev. Mod. Phys. 56, 255 (1984). 
[4] H. R. Miller and P. J. Witta, in Active Galactic Nuclei (Springer, Berhn, 1987), p. 202. 
[5] R. Fender, Annu. Rev. Astron. Astrophys. 42, 317 (2004). 

[6] P. Goldreich and W. H. Juhan, Astrophys. J. 157, 869 (1969); P. A. Sturrock, Astrophys. J. 

164, 529 (1971); V. L. Ginzburg, Sov. Phys. Usp. 14, 83 (1971); M. A. Ruderman and P. G. 

Sutherland, Astrophys. J. 196, 51 (1975). 
[7] F. C. Michel, Rev. Mod. Phys. 54, 1 (1982). 

[8] F. C. Michel, Theory of Neutron Star Magnetosphere (Chicago University Press, Chicago, 
1991). 

[9] V. S. Beskin, A. V. Gurevich, and Ya. N. Istomin, Physics of the Pulsar magnetosphere 

(Cambridge University Press, Cambridge, 1993). 
[10] M. Marklund and P. K. Shukla, Rev. Mod. Phys. 78, 591 (2006). 
[11] T. Piran, Phys. Rep. 314, 575 (1999); Rev. Mod. Phys. 76, 1143 (2004). 
[12] E. Tandberg-Hansen and A. G. Emshie, The Physics of Solar Flares (Cambridge University 

Press, Cambridge, 1988), p. 124. 
[13] M. L. Burns, in Positron-Electron Pairs in Astrophysics, eds. M. L. Burns, A. K. Harding, 

and R. Ramaty (AIP, New York, 1983). 
[14] V. I. Berezhiani, D. D. Tskhakaya and P. K. Shukla, Phys. Rev. A 46, 6608 (1992). 
[15] E. P. Liang et al, Phys. Rev. Lett. 81, 4887 (1998). 



9 



[16] C. Gahn et al, Appl. Phys. Lett. 77, 2662 (2000). 

[17] S. C. Wilks et al, Astrophys. Space Sci. 298, 347 (2005). 

[18] N. B. Narozhny et al., JETP Lett. 80, 382 (2004). 

[19] J. G. Lominadze et al. Phys. Scripta, 26, 455 (1982); M. E. Gedalin et al., Astrophys. Space 
Sci. 108, 393 (1985). 

[20] P. K. Shukla et al, Phys. Rep. 138, 1 (1986); M. Y. Yu, P. K. Shukla and L. Stenflo, Astrophys. 

J. 309, L63 (1986). 
[21] N. Iwamoto, Phys. Rev. E 47, 604 (1993). 

[22] G. Brodin, M. Markhmd B. EUasson and P.K. Shukla, Phys. Rev. Lett., 98 125001 (2007). 
[23] P. K. Shukla and L. Stenflo, Phys. Rev A, 30, 2110, (1984). 
[24] F. B. Rizzato, J. Plasma Phys. 40, 289 (1988). 

[25] V. L Berezhiani, L. N. Tsintsadze, and P. K. Shukla, J. Plasma Phys. 48, 139 (1992); Phys. 

Scr. 46, 35 (1992). 

[26] V. I. Berezhiani and S. M. Mahajan, Phys. Rev. Lett. 73, 1110 (1994). 
[27] V. L Berezhiani and S. M. Mahajan, Phys. Rev. E 52, 1968 (1995). 
[28] S. L Popel et al, Phys. Plasmas 2, 716 (1995). 

[29] P. K. Shukla, L. Stenflo and R. Fedele, Phys. Plasmas, 10, 310 (2003). 

[30] J. B. Kim et al, Phys. Lett. A 329, 464 (2004). 

[31] T. Cattaert, L Kourakis, and P. K. Shukla, Phys. Plasmas 12, 012319 (2005). 

[32] M. Marklund, P. K. Shukla, L. Stenflo, G. Brodin and M. Servin, Plasma Phys. Contr. Fus., 

47 25 (2005). 
[33] L. Stenflo, Phys. Scripta 14, 320 (1976). 
[34] P. K. Shukla, J. Plasma Phys. 72, 159 (2006). 
[35] G. Brodin and L. Stenflo, Phys. Scripta 37, 89 (1988). 

[36] L. Stenflo, P.K. Shukla and M. Y. Yu, Astrophys. Space Sci. 117, 303 (1985). 
[37] G. Brodin, Phys. Scripta, T113, 20 (2004); 

[38] V. A. Mironov, A. M. Sergeev, E. V. Vanin and G. Brodin, Phys. Rev. A. 42, 4862, (1990); 

G. Brodin, and J. Lundberg, Phys. Rev. E 57, 7041 (1998). 
[39] V. I. Karpman and H. Schamel, Phys. Plasmas 4, 120 (1997); V. I. Karpman, Phys. Scripta 

T75, 15 (1998). 



